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Truncated stellar disks in the near infrared. I. Observations 
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Abstract. We present NIR surface piiotometry of 11 edge-on galaxies obtained in the course of a long term project aimed 
at analysing the occurrence and type of the truncation of the outer disks. Observations were carried out at the 1.5 m CST 
(Carlos Sanchez Telescope) in Tenerife (Spain) using the CAIN infrared camera. 7 galaxies exhibit clear truncation on their 
disk profiles and 4 galaxies were observed to be clearly untruncated within observational limits. We describe the truncations as 
real, smooth and complete (as suggested by extrapolation and in the sense that the measured truncation curve goes into the noise 
at a truncation radius R„), following a decline proportional to (R„ - Ry" (where R is the radius). Despite its deep photometric 
reach, the data presented do not permit a detailed exploration of the region where optical data show a second slope. Special 
care was taken concerning the surface brightness deprojection of edge-on galaxies, which was carried out by two methods, one 
comprising the inversion of Abel's integral equation and the other following a numerical method. These methods gave nearly 
identical results. NIR observations of truncations could differ from observations in the optical, since the two domains trace 
different stellar populations. 
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1. Introduction 

It is widely accepted that the galactic disk brightness profiles, 
at least in the inner regions where the contribution of the bulge 
to the observed emission is still low, can be reproduced with 
an exponential function. At large radii, it is found that the stel- 
lar distribution often departs from this simple description and 
decreases much faster One of the goals of this paper is to con- 
tribute to determining the functional form of such a rapid de- 
cline. 

This decline of the stellar density was discovered by van der 
Kruit (1979), who denoted it as truncation. In this and subse- 
quent papers (van der Kruit & Searle, 1981a, 1981b) the basic 
properties of truncations were described as well as the first the- 
oretical ideas to explain them. Some recent studies have found 
that there is no complete truncation, but rather a sharp change in 
the exponential constant, called a "break" (Pohlen et al. 2002a, 
Pohlen et al. 2004 and references therein). 

Let us define the truncation curve as 



t(R) ^ fi(R) - fio(R) 



(1) 



as in Florido et al. (2001) (hereafter FOl). Here, R is the ra- 
dius, yu the actual surface brightness (in mag arcsec"^), and fio 
the same quantity obtained by the extrapolation of fi{R) in the 
inner disk, where the exponential decline can be followed un- 
ambiguously. 
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What is t(7?)? Does it depend on the observed wavelength? 
These questions, of fundamental importance to understanding 
this phenomenon, motivated the present paper and are only par- 
tially answered. Nevertheless, the observations reported here 
will contribute to achieving a final description, which is diffi- 
cult to assess because truncations (or "breaks") take place at 
a very low brightness in regions with a very low S/N ratio. 
The different theoretical models currently proposed are still far 
from obtain the functional form of t{R), but they should be able 
to do so eventually. Excellent reviews of the question have al- 
ready been made, such as those by van der ICruit (2001) and 
Pohlen et al. (2004). 

We opted to observe at NIR wavelengths, for the following 
reasons: firstly, the old stellar population, which makes up the 
great majority in the galactic disk of spiral galaxies, can be bet- 
ter traced, as it emits the bulk of its flux at NIR. There is also 
a contribution to the observed flux from the young population 
and from AGB stars, but this tends to be restricted spatially 
to specific areas within the disk (e.g. spiral arms and bulges) 
and does not change the basic trend in the stellar distribution 
throughout the disk. Another important argument is that the 
extinction effects are minimised and so the observed brightness 
distribution more closely resembles the true stellar distribution, 
while the stellar flux still largely dominates the measured flux. 
This is not the case at longer wavelengths, at which the ex- 
tinction is even lower. Also, the instrumental set-up allows us 
to reach outer regions on the galactic disks, where truncations 
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are normally found, due to a combination of an excellent dry 
site, a telescope designed and maintained for NIR observations 
and cryogenic instrumentation with background contamination 
controlled via a cold stop. 

Most previous studies (for example, Barteldrees & Dettmar 
1994; de Grijs et al. 2001; Kregel et al. 2002; Pohlen et al. 
2002a) have been carried out at optical wavelengths. To be 
complete, we mention the paper by Schwarzkopf & Dettmar 
(2000) as it contains a sample of 1 10 highly inclined galaxies, 
containing also K, data for 41 of the galaxies, but most of these 
are probably not deep enough. 

As a complement to this large amount of statistical data, 
detailed studies of the truncation curve are required. The ob- 
servations reported in the present paper and in FOl are con- 
siderably deeper than those of any previous study dealing with 
NIR photometry on the outskirts of spirals, here reaching val- 
ues of 20-23 mag arcsec"^ in J and 18.5-21 mag arcsec"^ in K^, 
which correspond to a face-on surface brightness in the range 
of 25-26 mag arcsec"^ in J and 24-25 mag arcsec"^ in K^. The 
increase of limiting magnitudes from edge-on to face-on views 
depends on the galaxy, mainly due to its inclination and thick- 
ness. These sensitivity limits enable us to extend the structural 
analysis to the outer areas of the disks, at galactocentric radii 
of comparable size to those achieved in optical studies. 

FOl presented NIR observation dealing with truncations, 
and the present paper can be considered as a continuation of 
the former. The telescope and the infrared camera proved to be 
very efficient for investigating the low surface brightness at the 
outskirts of disks. 

More NIR observations are necessary, especially taking 
into account that optical and NIR observations may differ sub- 
stantially and that this difference can be crucial. Recent star 
formation beyond the truncation (or break) radius would have 
a larger influence in the blue but very little influence in the NIR 
range. 

Different studies of truncations have used different tech- 
niques and so the results presented are not easily comparable. 
Some groups observe edge-on galaxies, while others consider 
face-on galaxies. Comparison is especially difficult when the 
data used are obtained at different wavelengths. For correct 
identification of the differences between optical and NIR obser- 
vations a parallel study of the same galaxies would be prefer- 
able, using the same analysis. This would be the best way to 
determine, definitively, wheffier there is a clear difference be- 
tween the results obtained using the different types of wave- 
lengths. We present here the results of our NIR survey which 
are noticeably deep and interesting in themselves. Observations 
at these wavelengths can provide vital clues to our understand- 
ing of the truncation phenomenon. 

2. The truncation curve in the literature 

Different descriptions for the truncation curve are classified in 
the schematic plot shown in Fig. 1. The face-on fi{R) and the 
truncation curve t(R) are plotted for four basic truncation types. 
This classification concerns the models considered to interpret 
the observed profiles. We have adopted this scheme of using the 
face-on fiir) in the truncation curve to homogenise data coming 
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Fig. 1. Different ways to describe the truncation phenomenon. 



from observations of both, face-on and edge-on galaxies. To 
render the last ones comparable with the former ones, we must 
deproject them. Deprojected profiles are intended to represent 
the true structure of the galaxy. 

Type 1. A sharp cut-off is introduced in the model. This ide- 
alisation is justified by the fact that the truncation curve interval 
is relatively narrow. 

The early observations by van der Kruit (1979) and van der 
Kruit & Searle (1981a,b) suggested profiles belonging to Type 
1 and were interpreted as a sharp cut-off in face-on galaxies that 
provided a relatively smooth truncation when the galaxy was 
seen edge-on. A sharp cut-off profile can appear as a smooth 
decreasing by the line of sight integration which is inherent 
to the observation of edge-on galaxies. Later on, Wainscoat et 
al. (1992) and many other works adopted this profile, after the 
work by van der Kruit and Searle. 

Type 2. A "smooth" profile, in which truncation is gradual, 
i.e. the first derivative is a continuous function and the trunca- 
tion becomes complete at Rtr, i.e. // = oo for /? > R,r (within 
observational errors). There is a transition region of a few kUo- 
parsecs. The edges are rounded in these profiles. 

Note that when we state that truncations are complete, 
this is an ideal description that is obtained by extrapolation. 
The zero intensity level is, of course, unreachable. However, 
in some galaxies the truncation curve is so well defined that 
"naked eye" extrapolation suggests that the truncation is com- 
plete. 

FOl presented NIR //(/?) and t{R) in agreement with this 
"smooth" description. The plots in their Fig. 3 for NGC 4013 
and NGC 5981 clearly suggest that the first derivative dfi/dR is 
a continuous function (in contrast, for example, to type 3). The 
real curve gradually moves away from the inner exponential to 
reach very large values. FOl proposed a function of the type 



TiR) 



constant 

iR,r-R)" 



(2) 



and after fitting the data, n was found to be close to unity. This 
function is a very simple one, a power law, for a continuous 
curve with a continuous derivative function. 

The type 2 profile is a smooth version of the classical 
sharp cut-off type 1 profile, avoiding a discontinuous derivative 
which is seldom found in nature. Profiles 1 and 2 are essen- 
tially the same, with the addition of a taring function, t(/?), to 
cormect the exponential disk and the sharp cut-off with slope 
00 in the type 1 profile. The introduction of a taring function 
has been considered by Casertano (1983) and by de Grijs et 
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al. (2001), who thus obtained a "soft cut-off" profile. However, 
they adopted other functions that differ from our equation 2. 

Type 3. The radial profiles show a two-slope shape sepa- 
rated by a rather sharp "break" radius. At the break radius, the 
inner exponential is followed by another exponential decline 
with a different slope. In this case, there is no real truncation, 
i.e. the truncation is not complete as in the type 2 profiles. This 
conclusion has also been reached by Pohlen et al. (2002a, 2004) 
and others. 

This profile is very well defined in NGC 5923 (Pohlen et 
al. 2004). It is noticeable that this galaxy is face-on (in early 
observations, truncations have been studied in edge-on galax- 
ies) which avoids deprojection effects. The two-slope profile 
description has also been reported by Erwin et al. (2005) and 
others. 

Type 4. The galaxy has no break and no truncation. This 
description corresponds to the early paper by Freeman (1970) 
where the disk was considered to be purely exponential. 

The profiles reported in this paper fit the type 2 "smooth" 
truncation profile, following Eq. 2. Our NIR data are not deep 
enough to permit a confident exploration, e.g. with sufficient 
S/N, of the region in which recent optical observations find the 
second slope after a break. 

In this paper, we present the data observed without theoret- 
ical preconceptions. It will be followed by a second one show- 
ing that the magnetic model matches the statistical properties 
derived from the data reported in this paper. 

3. Data 

As in FOl, the observations were carried out at the 1.5m CST 
at the Teide Observatory, Tenerife, with the NIR camera CAIN 
equipped with a 256^ NICMOS 3 detector array. The plate 
scale, with the wide field optics used, was 1.0"/pixel and the 
effective field of view was 4.3' x 4.3'. For the CAIN cam- 
era, the controller of the detector is based on a San Diego 
State University architecture, which provides low noise and 
high stability while permitting fast reading of the whole ar- 
ray, in about 50 milliseconds. In all cases, the read-out mode 
was the FOWLER mode, in which an equal number of non- 
destructive readings are taken, without time delays in between, 
immediately after the initial reset and before the final one, af- 
ter the selected integration time has elapsed. Thus, the pedestal 
and final flux levels are better determined, by averaging each 
group of readings independently. The final signal is then the 
difference between the last and first averages. The observa- 
tional characteristics are the same as in FOl (also explained 
in Castro-Rodriguez & Garzon 2003). The observations were 
carried out in 4 campaigns. The weather was generally good, 
being most of the nights photometric, on which we took the 
data from the standard stars and galaxies. A consistency check 
was performed for several bright objects against 2MASS point 
source catalogue, confirming that our measurements and those 
of 2MASS are consistent. 

The observed galaxies and their main physical parameters 
are depicted in Table 1 . Most of these were obtained from the 
LEDA database. The position angle (PA) is measured on the 
image with an automatic routine which delivers high precision 



values; for most of the galaxies the value obtained coincides 
with that of the LEDA database, the differences being < 2°. The 
orientation of the detector chip with respect to the sky was kept 
constant throughout the observation runs, so there was no sys- 
tematic effect due to changes in the instrument configuration. 
The PA was measured from the North celestial pole, increasing 
towards the East. Distances were obtained with Ho - 70 kms ' 
Mpc"^ taking into account the mean heliocentric radial velocity 
(km/s). All galaxies are edge-on, with an inclination / = 90° in 
LEDA. They were also selected to be, as far as possible, free of 
Milky Way stars in the outer edge and in such a way that their 
projected size {D25) fitted within the FOV. No other physical 
parameters (other than angular size and galactic latitude) were 
adopted for the selection; thus, although selection bias cannot 
be discarded, it is believed improbable. 

At these wavelengths, the night airglow emission is bright 
and rapidly varying; therefore, we used the classical observa- 
tional ON-OFF method with a nodding frequency compatible 
with that of the background variation in the sky and with an 
equal integration time on each pointing. The integration time 
interval was split into several frames to avoid exceeding the 
linear well depth of the detector due to the high sky flux, in 
particular in the Kj band. Alternating exposures were taken of 
the galaxy and the adjacent sky, separated by two minutes at 
most, for effective sky subtraction. Special care was taken with 
the flat-fielding and sky subtraction. For the reduction steps, 
we used several IRAF tasks and the IRAF package REDUCE, 
developed by R. Peletier. 

The observations were carried out in 4 campaigns. Table 2 
shows the references for each galaxy: Campaign 1 took place 
in October 2001, Campaign 2 in March 2002, Campaign 3 be- 
tween October and November 2002 and Campaign 4 in March 
2003. Table 2 also describes the following observational pa- 
rameters of the galaxies observed: the filter, the individual ex- 
posure time per frame (5 or 6 seconds for K, and H and 30 s for 
J), the number of frames at each position (4 for J and 20 or 24 
for Kj and H), the number of images per galaxy (usually 20 for 
J and 24 for K, and H); thus, the total exposure time was 40 (J) 
or 48 (K5 and H) minutes for the object. The last column shows 
the detection limit calculated as 3 times the STDDEV over the 
mean value of the background level, measured directly on the 
images in areas free from emission of the galaxy object. The 
typical seeing value was, on average, a little in excess of 1". In 
Fig. 4 we plot contour maps for these galaxies. 

4. Deprojection 

It is clear that the deprojection technique is critical in the reduc- 
tion process. In some cases, an initial inspection of the original 
profiles may suggest a truncation which disappears after depro- 
jection. To assure that the deprojection method works satisfac- 
torily, we carried out two tests: 

a) Deprojection was performed by two different, indepen- 
dent methods. The first of these is described in detail in FOl, 
while the second is based on Abel's integral equation (Binney 
& Tremaine 1987). Despite the fact that the numerical calcula- 
tions involved in these methods are intrinsically different, the 
results obtained by the FOl method and by using Abel's in- 
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Table 1. Basic properties of the galaxies. 



Galaxy 


RA 


DEC 


PA 


D 


Type 




log£)25 


log v,„ 




hh mm ss 


dd mm ss 


d 


Mpc 










(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


NGC 522 


01 24 45.91 


09 59 40.5 


33.3 


39.0 


Sbc 


-20.53 


1.44 


2.252 


NGC 684 


01 50 14.03 


27 38 44.2 


88.8 


50.4 


Sb 


-21.53 


1.53 


2.368 


MCG-01-05-047 


01 52 49.01 


-03 26 51.2 


161.4 


71.5 


Sc 


-21.77 


1.47 


2.410 


NGC 781 


02 00 09.02 


12 39 21.5 


13.0 


49.8 


Sab 


-20.87 


1.18 




NGC 2654 


08 49 11.91 


60 13 13.9 


65.0 


19.2 


SBab 


-20.09 


1.62 


2.295 


UGC 4906 


09 17 39.94 


52 59 34.3 


49.0 


32.6 


Sa 


-20.26 


1.30 


2.228 


NGC 2862 


09 24 55.10 


26 46 29.0 


114.0 


58.5 


SBbc 


-21.44 


1.41 


2.464 


NGC 3279 


10 34 42.61 


1111 50.7 


152.0 


19.9 


Scd 


-19.27 


1.44 


2.208 


NGC 3501 


11 02 47.35 


17 59 22.6 


28.0 


16.2 


Sc 


-19.05 


1.54 


2.147 


NGC 5981 


15 37 53.55 


59 23 30.9 


139.5 


36.1 


Sc 


-20.61 


1.43 


2.424 


NGC 6835 


19 54 33.09 


-12 34 02.5 


72.0 


23.0 


SBa 


-19.55 


1.38 


1.803 



Columns: (1) galaxy name; (2) and (3) coordinates (2000); (4) position angle; (5) distance, considering H„ = IQkms^^ Mpc^^; (6) 
morphological type; (7) absolute B-magnitude; (8) log of apparent diameter (D25 in 0.1') and (9) log of maximum velocity rotation (in km/s) 
from radio observations. All the parameters are obtained from the LEDA database 1 http://leda.univ-l yonl .fr^ except the position angles. 



Table 2. Observational parameters of the galaxies. 



Galaxy 


C 


B 


T„„/F 


F/P 


Pointing 


T 


3cr 








(s) 






ob.(m) 


mag/arcsec^ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 522 


3 


J 


30 


4 


20 


40 


20.6 




3 


K, 


5 


24 


24 


48 


21.1 


NGC 684 


1 


J 


30 


4 


20 


40 


20.5 




1 


K, 


5 


24 


24 


48 


18.6 


MCG-01-05-047 


3 


J 


30 


4 


20 


40 


22.7 




3 


K, 


5 


24 


24 


48 


20.7 


NGC 781 


3 


J 


30 


4 


20 


40 


22.9 




3 


K, 


5 


24 


24 


48 


20.7 


NGC 2654 


3 


J 


30 


4 


20 


40 


22.6 




3 


K, 


5 


24 


24 


48 


21.0 


UGC 4906 


4 


J 


30 


4 


20 


40 


22.0 




4 


Kj 


5 


24 


24 


48 


20.9 


NGC 2862 


3 


J 


30 


4 


20 


40 


22.7 




3 


K, 


5 


24 


24 


48 


20.7 




3 


H 


5 


24 


24 


48 


21.5 


NGC 3279 


2 


J 


30 


4 


20 


40 


22.0 




2 


K, 


6 


20 


24 


48 


20.4 


NGC 3501 


4 


J 


30 


4 


20 


40 


22.3 




2 


K, 


6 


20 


24 


48 


20.4 


NGC 5981 


2 


J 


30 


4 


20 


40 


21.9 




2 


K, 


6 


20 


24 


48 


20.2 


NGC 6835 


1 


J 


30 


4 


20 


40 


20.1 




1 


K, 


5 


24 


24 


48 


18.7 



Columns: (1) galaxy name; (2) campaign (1, oct. 2001; 2, march 2002; 3, oct.-nov. 2002; 4, march 2003); (3) band; (4) individual exposure 
time per frame; (5) number of frames at each position; (6) images per galaxy; (7) exposure time for the galaxy and (8) the 3cr level of the 
background variation over its mean value. 



tegration are almost indistinguishable. The FOl method was 
subsequently adopted, because of the few assumptions needed, 
and to use the same reduction process for the sake of coher- 
ence when merging the present data with those obtained in FOl . 
Thus, a larger and more homogeneous statistical sample is en- 
sured. Moreover, the FOl method can be generalised to more 
realistic conditions, including extinction. 



b) We have tested that the projection of the deprojected pro- 
file (x points) perfectly matched the original profile (A points) 
in a zero order test (Fig. 2). A deprojected profile corresponded 
to only one original profile and an original profile corresponded 
to only one deprojected profile. For example, our smooth de- 
projected profiles were not produced by a face-on profile with 
a sharp cut-off. 
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Fig. 2. Test for the deprojection method. Symbols are: A, input 
curve; □, deprojected curve; x, projection of the deprojected 
curve. 



We took the z-integrated profiles as the original input pro- 
file for the deprojection method, i.e. those profiles obtained 
by summing all the input radial profiles for each z. These z- 
integrated profiles were obtained as follows: first the galaxies 
were rotated to make the major axis perpendicular to the or- 
dinate axis. From I(r,z), where I is surface brightness, r is the 
abscissa and z the ordinate, the z-integrated profile is then ob- 
tained as /(r) - lll^Jir, z), where a and b are the limits of a 
box which contains the 3cr isophote of galaxies. We have not 
masked any foreground stars. We thus create an artificial pro- 
file. If we wished to obtain the truncation curve as a function of 
z, this z-integration would prevent it. This could be interesting, 
for instance, to see if the truncation were difl'erent for the thin 
and the thick disks, which is especially important in early type 
spirals. Moreover, in early type galaxies the extended bulge 
could smear out a possible truncation in a z-integrated profile. 
This cannot be undertaken when using face-on observations. 

Certainly, by z-integrating we cannot obtain the emission 
per unit volume in the galaxy. However, this procedure has 
two important advantages: firstly, we can work with a profile 
which has a much better signal/noise ratio, and secondly, the ef- 
fect of the large extinction dependence on z is greatly reduced. 
Suppose, for instance, a dust lane for z=0 or around z—0. The 
row z=0 would give meaningless results because of severe ex- 
tinction. By summing all the z rows, the effects of localized 
dust lanes are reduced. 

The z-integration apparently leads to a loss of information. 
However, this efl'ect is irrelevant if the final product of the de- 
projection method is the face-on view. In Fig. 3 we plot both the 
z-integrated and the deprojected profiles for a selected galaxy, 
NGC 2862, for the three filters to illustrate the difl'erence. 
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Fig. 3. NGC 2862: z-integrated profile (upper) and deprojected 
profile (lower) for each filter. 



It might be thought that in order to obtain the face-on view 
the sequence "deprojection, then z-integration" is not equiva- 
lent to the sequence "z-integrating, then deprojection". From 
the description of the deprojection method in FOl it can be 
demonstrated that the two sequences are equivalent. This is 
shown in Appendix A. 

Therefore, the method itself does not introduce impor- 
tant errors (as concluded in FOl). However, important errors 
could be introduced if the required assumptions are not accom- 
plished. These assumptions are axisymmetry, the absence of 
extinction, and an inclination equal to 90". In other words, the 
deprojected profile need not necessarily coincide with the real 
face-on profile. They only coincide when the required assump- 
tions are met. 

With respect to axisymmetry, it is clear that when the 
galaxy lacks this property then deprojection cannot reproduce 
the real face-on view. This source of potential errors is in- 
evitable in this and in any other deprojection method. The 
problem has infinite solutions, as from a 2D distribution it is 
not possible to obtain a 3D one without certain assumptions. 
Fortunately, in the galaxies observed, the noticeable east-west 
symmetry suggests that they are reasonably axisymmetric. 

An estimation of this error is not easy. NIR data of face- 
on galaxies have not been reported in the literature up to these 
large radii, so we cannot use additional real data for this es- 
timate. The only galaxy for which we have at our disposal 
NIR data at large radii is the Milky Way. Wainscoat et al. 
(1992) have modelled the Milky Way in order to reproduce 
IRAS data, including NIR arms which are the mean sources 
of non-axisymmetries. The over-density due to one NIR arm 
is very small. From these data about the relative importance 
of arm with respect the unperturbed disc and taking into ac- 
count the contribution to the flux of each spectral type, from 
Porcel (1997) we obtain that the contribution is mainly due to 
stars of spectral types B8 V to F5 V and is only 1%. When the 
propagation of errors in the deprojection method is considered, 
this figure becomes representative of the error introduced by 
one-arm asymmetries by the deprojection method. Galaxies re- 
ported here could be less symmetric than the Milky Way, but 
we estimate that this kind of error can be ignored. 

Extinction is another inevitable source of errors, partic- 
ularly if there is a non-axisymmetric distribution of dust. 
However, the problems introduced by extinction are reduced 
in the present approach because we work in the NIR and be- 
cause prior z-integration avoids the efl'ects of inhomogeneities 
on the distribution of dust, e.g. dust lanes in the plane. 
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We used only galaxies with 1=90", following LEDA. 
However, the LEDA galaxies may deviate typically 2° from 
edge-on, which could be a source of errors. Our galaxies have 
a minor-axis of 70" sin 2° = 2.4", which is higher than our 
resolution about 1". In Appendix B we show that inclinations 
i > 88° do not introduce noticeable errors. 

Another important step in the analysis was to calculate 
the slope of the exponential function of the form I{R) = 
Iq exp(-R/h) in the inner part. This parameter h is present in 
the definition of TiR) and it is necessary to be able to estimate 
Rtr/h, another parameter that is derived in all the studies made 
of truncations/breaks. The region for calculating h was decided 
upon by visual estimation, h is calculated in the deprojected 
profiles. Noise and intrinsic variations in surface brightness do 
not permit the determination of h with an error of less than 
about 10%. 

The use of a visual decision regarding the interval for cal- 
culating h does not introduce large uncertainties, because of 
the broad region in which the exponential behaviour provides a 
good fit. If the photometric profiles are divided into bulge, ex- 
ponential and truncation shapes (with no sharp transitions) the 
choice of an excessively large interval would provide an under- 
estimation of h, while the choice of one that was excessively 
small would give large errors. But excluding these two extreme 
cases, the value of h is nearly independent of the choice of the 
interval. Therefore, this subjective choice does not introduce 
additional errors. 

The main results are presented in Table 3, for each side 
of the galaxies. After the deprojection process, by fitting the 
exponential profile, we obtain jjg, the central surface brightness 
of the disk, and h, the radial scale length. Then we estimate Rtr, 
the truncation radius, and n, the exponent in equation (2). We 
also calculate Rtr/h. 

5. Results 

The galaxies observed can be grouped taking into account the 
existence of truncation, as follows: 

a) Truncated galaxies: NGC 522, MCG-01-05-047, NGC 
781, NGC 2862, NGC 3279, NGC 3501 and NGC 5981. Our 
study is mainly focused on this group. 

b) Untruncated galaxies: NGC 684, UGC 4906, NGC 2654 
and NGC 6835. Of course we cannot assure that these are not 
truncated, but only state that no truncation was observed for 
face-on values of J<23 mag/arcsec^ in NGC 684, for J<24 
mag/arcsec^ in UGC 4906, for J<25 mag/arcsec^ in NGC 2654 
or for J<22 mag/arcsec^ in NGC 6835. We refer to the J fil- 
ter limits since in these galaxies the J band data covers a more 
extended region in the object, as can be seen in Fig. 4. In con- 
sequence we did not detect truncations out to 3.5 radial scale- 
lengths for NGC 684, 5.0 for UGC 4906, 5.9 for NGC 2654 or 
3.3 for NGC 6835. 

We deduce that approximately 2/3 of the considered galax- 
ies are truncated, a sufficiently large value to highlight, once 
more, the importance of this effect. Although this figure must 
be approached with some caution, the overall conclusion must 
be that truncations/breaks are a fairly common feature of nor- 
mal spirals. It is also evident that we cannot speak of untrun- 



cated galaxies, but merely of galaxies for which truncation is 
not observed within the limits of our observational capabilities. 

In Fig. 5, we plot the deprojected profiles for these galaxies. 
Fig. 5 also show t(R) obtained from Eq. 1 for these galaxies and 
for two galaxies (MCG-01-05-047 and NGC 5981), included as 
examples, we show the fitted truncation curve. 

Table 3 shows that the truncation radii obtained with the J 
filter is always larger than the value obtained with the K, filter, 
the ratio being RtriJ)/RtriKs) = 1.15. The dispersion of the 
data is 0.04, but our error in detemnining this quantity is larger, 
at about 0.5. 

The identification performed is not free of definition prob- 
lems. Let us assume a galaxy with an extended (not sharp) 
"break" region and with the second outer exponential lying be- 
yond observation capabilities. This could have been classified 
as "smooth" despite being a "break" galaxy. The data for some 
of our galaxies strongly suggest that the "smooth" description 
is physically different from the two-slope type, but this identifi- 
cation is not easy. Subjective inspection is probably essential to 
determine whether the change in the slope is gradual or sharp 
or, equivalently, whether the first derivative is a continuous or a 
discontinuous function. Given the mounting body of evidence 
provided by Pohlen et al. (2002a, 2004), Perez (2004), Erwin 
et al. (2005), Trujillo & Pohlen (2005) and others in obtain- 
ing clear two-slope profiles, we could have fit the data to this 
model. But considering that our data do not permit the explo- 
ration of the second-slope region and considering that the NIR 
and the optical profiles do not necessarily coincide, we follow 
an independent approach. 

In the previous paper by FOl, the deprojected profiles of 
NGC 4013 (both J and K„ both sides), NGC 4217 (K„ left side) 
and NGC 5981 (K,, right side) are clear examples of a gradual 
decline and a continuous first derivative. Notice that NGC 5981 
is present in both papers (in the first one only in K,). The K,- 
images are very similar as expected, but in the first campaign 
better observing conditions allowed us to reach slightly more 
distant radii. The difficulties inherent in subjective discrimina- 
tion cannot be easily surmounted by statistical methods, as the 
break point should be determined manually. Thus, we prefer 
the "smooth" truncation, type-2 profiles with gradual decline 
and gradual variation of the first derivative. 

We then fitted the observation points to the function in Eq. 
2. Rtr was obtained by extrapolation and the "constant" and the 
value of n by standard fitting methods. The n value is particu- 
larly interesting. The results are shown in Table 3. 

This fitting features the compromise of a gradual transition 
from a pure exponential to a continuously bending profile. This 
is, however, provided automatically by the function given in 
Eq. (2). TiR) becomes negligible for values of R not far from 
Rtr- For example, assuming n=l as a typical value, it is found 
that T is only 10% of the value for R = 0.95R,r, for a value of 
R one radial length less than Rtr. Therefore, the result is prac- 
tically independent of the starting point in the fitting (provided 
the bulge zone is excluded). The real curves in Fig. 5 show how 
fast T decreases far from the truncation radius. 
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Table 3. Photometrical and structural parameters of the target list 



Galaxy 


Passband 


oiue 




rl 


n 






n 










(mag/arcsec ) 


(kpc) 


(arcsec) 


(kpc) 


(arcsec) 








(^) 


(3) 


(4j 


(->) 


(o) 


(') 


(s) 


(9) 


(10) 


NGC 522 


J 


SW/NE 


20.6/20.6 


7.2/6.1 


38/32 


14/14 


73/72 


1.0/1.0 


1.9/2.2 




K, 


SW/NE 


21.1/20.9 


8.0/6.3 


42/33 


11/11 


59/60 


0.5/0.7 


1.4/1.8 


NGC 684 


J 


NW/SE 


18.8/18.6 


5.3/4.6 


22/19 














NW/SE 


16.8/17.3 


3.7/4.1 


15/17 










MCG-01-05-047 


J 


NW/SE 


20.9/20.9 


10.4/11.3 


30/33 


30/29 


85/84 


0.6/0.7 


2.8/2.6 




K, 


NW/SE 


19.5/19.6 


8.4/9.4 


24/27 


28/27 


82/79 


1.0/1.0 


3.4/2.9 


NGC 781 


J 


SW/NE 


19.4/19.8 


2.6/3.1 


11/13 


11/11 


46/45 


0.9/0.8 


4.2/3.5 




K. 


SW/NE 


18.3/18.0 


2.6/2.3 


11/9 


10/10 


41/40 


0.5/0.3 


3.5/4.3 


NGC 2654 


J 


NE 


18.9 


1.5 


17 












K, 


NE 


17.3 


1.4 


15 










UGC 4906 


J 


SW/NE 


18.3/18.3 


1.9/1.7 


12/11 












K, 


SW/NE 


17.3/17.3 


1.7/1.7 


11/10 










NGC 2862 


J 


NW/SE 


19.5/19.6 


5.9/6.0 


21/21 


22/22 


iipe 


1.1/0.9 


3.7/3.6 




K, 


NW/SE 


18.5/18.5 


5.9/5.6 


21/20 


19/18 


66/65 


0.8/0.7 


3.2/3.3 




H 


NW/SE 


18.7/18.7 


6.4/5.8 


23/20 


21/20 


74/69 


1.1/0.3 


3.3/3.4 


NGC 3279 


J 


NW/SE 


19.9/19.6 


4.7/3.1 


48/32 


8/8 


79/82 


1.5/1.0 


1.6/2.6 




K, 


NW/SE 


18.5/18.3 


3.8/2.9 


39/30 


7/7 


69/74 


1.2/1.2 


1.8/2.4 


NGC 3501 


J 


SW/NE 


20.2/19.9 


2.1/1.7 


27/22 


8/9 


100/110 




3.7/5.1 




K. 


SW/NE 


18.1/18.0 


1.6/1.6 


21/20 


6/ 


82/ 


0.5/ 


3.9/ 


NGC 5981 


J 


NW/SE 


19.4/19.3 


4.3/4.2 


25/24 


15/14 


88/79 


0.9/0.9 


3.6/3.3 




K. 


NW/SE 


18.1/18.0 


4.0/3.7 


23/21 


13/12 


74/70 


0.9/0.0 


3.2/3.3 


NGC 6835 


J 


NW/SE 


18.1/18.1 


1.8/1.8 


16/16 












K, 


NW/SE 


17.4/17.4 


1.8/1.8 


16/16 











Columns: (1) galaxy name; (2) filter; (3) side of the galaxy; (4) central surface brightness; (5) & (6) radial scale-length in kpc and in arcsec; 
(7) & (8) truncation radius in kpc and arcsec; (9) parameter n obtaining from Eq. (2); (10) ratio between truncation radius and scale-length. All 
values are deprojected (face-on). 



The mean value of n is 0.8 with cr - 0.3. This supports the 
conclusion by FOl who also found a coefficient close to unity. 
Therefore, we propose 



Notice that the face-on and edge-on truncation curves differ 
but Rtr is the same (Fig. 3). Therefore, this value is unalfected 
by the use of a deprojection technique. The Rtr values obtained 
in other studies without deprojection are equally valid for sta- 
tistical studies. 

The Rtr and n parameters are obtained independently. The 
value of R,r tells us "where" the truncation takes place and the 
values of n tell us "how sharp" the truncation is. 

6. Conclusions 

In order to obtain more information about the truncation of 
the stellar disk, it is necessary to incorporate observations in 
the NIR, which is the best tracer of the old population and 
relatively free of extinction problems, and to use larger tele- 
scopes. On the other hand, observations in the optical and in the 
NIR should be considered complementary, as both wavelength 
ranges may give information about different aspects of the fea- 
ture. In fact, they may differ not only in the value of certain 
fitting parameters, but in the functional form of the truncation 
curve itself. For example, beyond the truncation (or break) ra- 



dius, non negligible star formation could affect the young pop- 
ulation distribution but not that of the old population. If stars 
did not move from their birthplace, a difference in the blue and 
the NIR profiles would not be expected, but stars can move. 

If there is a (smooth and complete) truncation in the IR, as 
shown by FOl and this study, and a sharp break in the optical, 
as shown by Pohlen et al. (2004), this is not necessarily a con- 
tradiction; it could be the sign of differing behaviour patterns 
in the optical and in the near infrared, for dynamic reasons. We 
choose the traditional description introduced by van der Kruit 
(1987), introducing a taring function to smooth out the sharp 
cut-off, also called a soft cut-ofF, and propose a truncation curve 
proportional to (Rir - RY", with n - 0.8 ± 0.3, a value close 
to 1, i.e. it becomes complete alR - Rtr- This is the profile of 
type 2, labelled "smooth" in Fig. 1. 

The type 3 profile is clearly different from that of type 2, 
either due to the different information provided by the NIR and 
the optical, or due to insufficient coverage of the NIR observa- 
tions. In the light of the former consideration, we have chosen 
the simplest fit suggested by our data, but they are not deep 
enough to conclude that an external region with a second slope 
does not exists. 

Unfortunately, there are not many galaxies measured in 
both the optical and the NIR. One of them is NGC 522, which 
was measured in the optical by Pohlen et al. (2002b). This 
galaxy is a very good example of a two-slope with a nar- 
row break region in between, when measured in the optical. 
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However, as can be seen in Fig. 5, the J-profile on the right hand 
side seems to be smooth. The left side as well as the right Ks 
side could be described as a two-slope. The data does not allow 
us to clearly distinguish between these two types. Hence, no 
clear conclusions can be obtained from this comparison. This 
important point should be addressed in a future study. 

As the radial scale height is a natural unit of length in a spi- 
ral galaxy, Rtrlh is a parameter characterising truncations. For 
these galaxies, we found the following values for Rtrlh: 3.2 for 
J and 2.9 for K„ with cr - 0.9 and 0.8 respectively for each 
filter. The colour dependences of Rfr will be discussed in Paper 
II. These values are in reasonable agreement with previous es- 
timations. For comparison, it should be taken into account that 
the diff"erence between the radius at which the truncation curve 
leaves the inner exponential and the radius at which the trun- 
cation finishes is of the order of 20-30", and that in the break 
description only the first value is meaningful. 

It is not easy to establish how often galaxies are truncated, 
due to the difliculty in determining whether a given galaxy is 
non truncated or whether Rtr is too distant to be observed. If 
we consider 7 out of 11 galaxies truncated, this is equivalent 
to approximately 2/3. However, this figure should be consid- 
ered with some caution as our sample is small. However, we 
confirm that the frequency of truncated galaxies is high, which 
makes this topic of crucial importance in our understanding of 
the evolution of disk galaxies. 
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Fig. 7. High inclination view of the galaxy. 



Fig. 6. Definition of the involved quantities in Appendix A. 



Appendix A: z-integrating + deprojecting = deproject- 
ing + z-integrating. 

For our purpose the two reduction steps of z-integrated and 
deprojecting can be interchanged, as we are only interested in 
obtaining the face-on view of the galaxy. 

We observe I{r,z), the surface brightness of the edge-on 
galaxy. We want to obtain /(/?), the surface brightness as a func- 
tion of the galactocentric radius, for a face-on galaxy. 

We take the pixel size as unity. 

We can follow two procedures: 

a) First summing in z, then deprojecting 

b) First deprojecting, then summing in z. 

In the first case, from I{r, z) we obtain I{r) = "LJir, z), then 
obtain I{R) - j^-;r^ (quantities are those defined in FOl). This 
last equation is schematic. Of course, I{R) is obtained by a se- 



ries of more complicated steps but they are essentially of this 
type and to foUow the series step by step does not introduce 
any change. 

In the second case from /(r, z) we obtain /(/?, z) as /(/?, z) = 
(observe that the areas A depend on R and r but not on z). 
Then we sum in z 



/(/?)= S,/(/?,z)=E, 



1 



I(r,z) ^ 

AiR, r) A{R, r) 



^Ar,z) = 



1 



A{R, r) 



obtaining exactly the same result as in the first procedure. 

Appendix B: An inclination 90" >i>88'' does not intro- 
duce important errors. 

These errors depend not only on i but also on the angular 
resolution. Suppose a galaxy with / = 88°. Along the minor 
axis distances are reduced by a factor sin 2" - 0.035. We want 
to know the size OO' perpendicular to the major axis (Fig. 7). 
In an exponential disk, the galaxy has no end and 00' - oo. 
However we want a characteristic value as at P the density is 
higher and decreases rapidly when we go from P to O'. We 
must find I{x) for a given r. 

We assume / - IoCxp(-R/h), hence /(x) = 
Ioexp-{r^ + x^y^^, where r is a constant. First we work 
with unprojected distances. By a Taylor expansion 



1/2 



if! 
2~r 



therefore 




I{r) Fig. 8. Definition of the involved quantities in Appendix B . 



istic OO' distance, considering that it corresponds to a value of 
X at which I{x) is I{P)le 

_i£ 1 

g 2 rh = — 

e 

2rh 

x= ^/2rh 



For typical values in our galaxies, r = 40", h = 10", we 
have X = 30". And now we project 30" sin 2° - I". This is of 
the order of the pixel size. Therefore our pixel size is too large 
to appreciate the i + 90° inclination. Therefore, variations in z 
are mainly attributable to real variations in the vertical direction 
for i > 88". 



where I{P) - I(R = r) is a constant. We now find a character- 
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Fig. 4. Contour maps for the observed galaxies for J, H and Ks filters. The isophotes are equidistant (in units of n*3cr equiv. to a 
step of +0.75 mag/arcsec^) starting at our noise level (see Table 2). Galaxies are rotated, but N is closer to the top and E to the 
left. 
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Fig. 5. Deprqjected profiles for all galaxies. Left: The solid line is the exponential fit with the parameters in Table 3. Horizontal 
bars indicate the region in which the radial scale length was calculated and the arrow the truncation radius. (T) is for truncated 
galaxies and (U) for untruncated ones. Right: For each of these galaxies we have represented the truncation curve (r) from the 
equation 1 and, as an example, we have plot the fitted function of Eq. 2 for two of them. 
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Fig. 4. continued 
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